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We study the Landau gauge gluon propagator D{p) in the 3d SU{2) lattice gauge theory. We 
show that in the infinite-volume limit the expectation values over the Gribov region Q are different 
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conclusion does not change when spacing a tends to zero. 
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I. INTRODUCTION 

There are various scenarios of confinement based on 
infrared behavior of the gauge dependent propagators. 
For example, in the Gribov-Zwanziger (GZ) confine- 
ment scenario [l|, I4] the Landau-gauge gluon propa- 
gator D{p) at infinite volume is expected to vanish in 
the infrared (IR) limit p ^- 0. At the same time, a re- 
fined Gribov-Zwanziger (RGZ) scenario '3:^-5] allovi^s a 
finite nonzero value of D(0). The nonperturbative lat- 
tice calculations are necessary to check the validity of 
each scenario as well as to check the results obtained 
by analytical methods, e.g., the (truncated) Dyson- 
Schwinger equations (DSE) approach. The DSE scal- 
ing solution predicts that the propagator tends to zero 
in the zero-momentum limit [y, |7[ in accordance with 
the GZ-scenario. Another - decoupling - solution [8- 
lll| allows a finite nonzero value of £'(0) in conformity 
with RGZ-scenario. 

The 3d SU{2) theory can serve as a useful test- 
ground to verify these predictions. It is also of in- 
terest for the studies of the high-temperature limit of 
the 4d theory. Last years the 3d theory has been nu- 
merically studied in a number of papers |12hl8]. It 
has been shown that the propagator has a maximum 
at momenta about 350 -^ 400 MeV and that zero mo- 
mentum propagator -D(O) does not tend to zero in the 
infinite- volume limit. 

The Gribov copy problem still remains one of 
the main difficulties in computation of the gauge- 
dependent objects (for 3d case see, e.g., [l^ and ref- 
erences therein). 

The manifold consisting of Gribov copies provid- 
ing local maxima of the gauge fixing functional F{U) 
(defined in Section |TT1) and a semi-positive Faddeev- 



Popov operator is termed the Gribov region Q, while 
that of the global maxima is termed the fundamental 
modular region (FMR) F C r2 [l^ . Our gauge- fixing 
procedure is aimed to approach F. 

In paper [20| it has been claimed that although 
there are Gribov copies inside Gribov region f2 , they 
have no influence on expectation-values in the ther- 
modynamic limit i.e. for any gauge noninvariant ob- 
servable O 



{0)n = {0)i 



(1) 



In our recent paper 'iSJ we attempted to check this 
statement. We calculated gluon propagators D{p) on 
different lattices (for p = as well as for p ^ 0) and 
then extrapolated the values of D in the thermody- 
namic limit. It has been shown that in the thermo- 
dynamic limit L — > 00 the value of propagator -D(O) 
clearly depends on the choice of the gauge copy. 

The most of our calculations in [1^ has been per- 
formed at /3 = 4.24 (a = 0.168 fm). The main goals of 
this paper are (a) to find confirmation of our observa- 
tions made in [l8| employing different (larger) values 
of (3 and (b) to draw some definite conclusions about 
the continuum limit of the theory. 

In Section In] we introduce the quantities to be com- 
puted and give some details of our simulations. In 
Section IlIII we present our numerical results. Conclu- 
sions are drawn in Section [IVI 



II. 



MAIN DEFINITIONS AND DETAILS OF 
THE SIMULATION 



We consider 3d cubic lattice L^ with spacing a. To 
generate Monte Carlo ensembles of thermalized con- 



figurations we use tlie standard Wilson action 



S^(3 



Xjfj^yv 



1- 



2 ^"^\y x^iU x+jla:uU ^j^0g^.^^U ^^ 



(2) 

where j3 — A/g^a, /i is a vector of unit length, along 
the /zth coordinate axis and gs denotes dimensionful 
bare coupling. [/^.^ G SU{2) are the link variables 
which transform under local gauge transformations g^ 
as follows: 



^<°' - Iv^) 



5^5^aj(o)4(o)). (8) 



tj.=i &=i 



In what follows we use the gluon propagator D{p) 
normalized at fi = 2.5 GeV, so that p^D{p) — 1 for 



2 2 

P = M • 



We employ the usual choice of the Landau gauge 
condition on the lattice 1221 



U.f. 4 U^^ = glUx^9.+f.a , g. e SU{2) . (3) {dA)ix) ^ i J^ U.ix + ^) - A,ix ^ ^)) ^ 



In Table |I] we provide the full information about the 
field ensembles used in this investigation. The scale 
is set in accordance with |2l| where string tension is 
V^ = 440 MeV. 

We study the gluon propagator 



Dtlil) 



L3 



'^e:cp (^iqx + — qifi - i>)j (4) 



{Alix + y+^)A:{y+^)), 



where the vector potentials A'f^{x) are defined as fol- 
lows in : 



A, 






t<-^-i-,iU.^.-UlJ^ (5) 
6=1 ^ 9Ba 



and the momenta q^ take the values g^ = 27rn^/aL, 
where n^ runs over integers in the range — L/2 <n^< 
L/2. The gluon propagator can be represented in the 
form 



D^M) 



S'^Sf^.DiO), 



2 „;„ 9^0 _^ „2 _ Y-3 9 






where Pp. = — sin -^ and p^ = J2n=i pfi- For p 7^ 



a 
one arrives at 



^(P) 



1 1 



6 (La) 



^^(i^(g)i^(-g)), (6) 



M=l ft=l 



where 



i^(g) = a^^4(x + ^)exp(zg(. + ^)), (7) 



and the zero-momentum propagator has the form 



a ^^ \ ■ 2 ' '^' 2 

(9) 
which is equivalent to finding a local extremum of the 
gauge fixing functional 



Fuig) 



L y^ i Tr [/» 
7-3 Z^ 9 ^x^l 



3L3 



(10) 



Xfl 



with respect to gauge transformations g^^. 



To fix the gauge we choose for every gauge orbit a 
representative from T |19| , i.e. the absolute maximum 
of the gauge fixing functional F{U). This choice is well 
consistent with a non-perturbative PJLZ gauge fixing 
approach [2J, [2J| which presumes that a unique repre- 
sentative of the gauge orbit needs the global extremum 
of the chosen gauge fixing functional. Also in the case 
of pure gauge C/(l) theory in the Coulomb phase some 
of the gauge copies produce a photon propagator with 
a decay behavior inconsistent with the expected zero 
mass behavior |25l - l27| . However, the choice of the 
global extremum permits to obtain the massless pho- 
ton propagator. 

For practical purposes, it is sufficient to approach 
the global maximum close enough so that the system- 
atic errors due to nonideal gauge fixing (because of, 
e.g., Gribov copy effects) are of the same magnitude 
as statistical errors. This strategy has been checked 
in a number of pap ers on Ad and 3d theory studies 
for both SU{2) 51 Hi-Ill and 5'f/(3) [H ^ gauge 
groups. 

The gluon propagator in the deep infrared region 
can be reliably evaluated only when the effects of 
Gribov copies are properly taken into account. The 
gauge-fixing algorithm which we use was already suc- 
cessfully employed in the Ad theory at both zero 
|30l . l3l| and nonzero 32, ^5s] temperature. There are 
three main ingredients of this algorithm: powerful 
simulated annealing algorithm, which proved to be ef- 
ficient in solving various optimization problems; the 
flip transformation of gauge fields, which was used to 
decrease both the Gribov-copy and finite-volume ef- 
fects |30l - l32| : simulation of a large number of gauge 



P = 7.09 (a = 0.094 fm) 
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P = 10.21 (a = 0.063 fm) 
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TABLE I: Values of lattice size, L, number of measure- 
ments rimeaa and number of gauge copies Ucopy used 
throughout this paper. 



copies for each flip sector in order to further decrease 
the effects of Gribov copies. 

All details of our gauge fixing procedure can be 
found, e.g., in ]\^. For readers convenience, we will 
describe it shortly here. 

Firstly, we extend the gauge group by the transfor- 
mations (also referred to as Z2 flips) defined as follows: 



U{u.,^) 



- u. 



x,^ 



if 11 = 1/ and 



Ux^^ otherwise 



which are the generators of the Z| group leaving the 
action ([2]) invariant. Such flips are equivalent to non- 
periodic gauge transformations. A Polyakov loop di- 
rected along the transformed links and averaged over 
the 2-dimensional plane changes its sign. Therefore, 
the flip operations combine the 2^^ distinct gauge or- 
bits (or Polyakov loop sectors) of strictly periodic 
gauge transformations into one larger gauge orbit. 

We use the simulated annealing (SA), which has 
been found computationally more efficient than the 
use of the standard overrelaxation (OR) only [23, [33, 
ISq . The SA algorithm generates gauge transforma- 
tions g{x) by MC iterations with a statistical weight 
proportional to eyi\i{iV Fij[g\/T) . The "tempera- 
ture" T is an auxiliary parameter which is gradually 
decreased in order to maximize the gauge functional 
_?[/[(/] . In the beginning, T has to be chosen suf- 
ficiently large in order to allow traversing the config- 
uration space of g{x) fields in large steps. T is 
decreased with equal step size. The final SA temper- 
ature is fixed such that during the consecutively ap- 
plied OR algorithm the violation of the transversality 



^^ max I (a^^) (a;) I < e 

2 X,C 



(11) 



decreases in a more or less monotonous manner for 
the majority of gauge fixing trials until the condition 
(fTTj) becomes satisfied with e = 10~^. 

To finalize the gauge fixing procedure we apply the 
OR algorithm with the standard Los Alamos type 
overrelaxation. In what follows, this method is la- 
belled FSA ("Flipped Simulated Annealing"). 

We then take the best copy {he ) out of many gauge 
fixed copies obtained for the given gauge field config- 
uration, i.e., a copy with the maximal value of the 
lattice gauge fixing functional Fjj as a best estimator 
of the global extremum of this functional. 

To demonstrate the effect of Gribov copies, we also 
consider the gauge obtained by a random choice of 
a copy within the first Gribov horizon (labelled as 
"/c " — first copy), i.e., we take the first copy obtained 
by the FSA or SA algorithms. It is instructive also 
to compare be and /c propagators with the worst copy 
{wc ) propagators which correspond to the choice of 
the gauge copy with minimal value of the gauge fixing 
functional. 



III. NUMERICAL RESULTS 

To estimate the infinite- volume limit L — > 00 of the 
zero-momenturn gluon propagator D{Q\ L) we apply 
the fit-formula |14| 



D(0;L)==ci + C2/L 



(12) 



In Fig. [T] we show our values of 15(0; L) calculated 
for /3 = 7.09 (left panel) and /3 = 10.21 (right panel) 
for be , /c and we calculated by FSA method. Straight 
lines represent fits according to Eq. (IT2t . 

In all cases the values of D(0) in the thermodynamic 
limit differ from zero which is in agreement with the 
statement made in [ij, |l^ . 

Moreover, Fig. [T] demonstrates another interesting 
phenomenon : the Gribov copy infiuence remains 
rather strong even in the thermodynamic limit. In- 
deed, the infinite-volume extrapolation of D^'^ifi) dif- 
fers from infinite- volume extrapolation of D^'^{0) (as 
well as for £'""^(0)) which agrees with our observation 
made in [3 for /? = 4.24. 

Therefore, the expectation values over the Gribov 
region 51 are different in the infrared from that calcu- 
lated in the fundamental modular region F that dis- 
agrees with the statements made in |2(|] . This is main 
result of our paper. 

As in [H] we found that while for finite L D)?^{Q, L) 
is higher than D^psj^{^\ L), in the infinite- volume they 
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FIG. 1: D(0) as a function of 1/aL for fi = 7.09 (Left) and for ^ = 10.21 (Right). 
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FIG. 2: Vy(0) as a function of 1/aL for /3 = 7.09 (Left) and fi = 10.21 (Right). 



coincide within errorbars. This remarkable agreement 
confirms the rehability of our estimation of D(0) in the 
infinite- volume limit. 

For better illustration of our main result we cal- 
culated additionally the averaged difference between 
/c and he propagators normalized to D^'^{p;L = oo) 



Wip) 



Pf^jp) - D'"'{p) 
D'"'{p]L = oo) 



(13) 



In Fig. [2] we show the dependence of 14^(0) on the 
inverse lattice size both for FSA and SA algorithms. 
For FSA the value of 14^(0) decreases for rising size, 
while for SA W{0) grows. 

For smaller volumes the values of W{0) for SA al- 
gorithm, Wsa{Q), are close to zero ("Gribov noise") 
while Wfsa{0) is at its maximum. The last corre- 
sponds to strong effects of flip-sectors (see, e.g., |18jV 
However, with increasing volume Wsa{0) is increas- 
ing indicating the increasing role of the copies within 
given flip-sector. In opposite, decreasing of Wfsa{0) 



with increasing volume implies that the role of the 
flip-sectors reduces. 

Remarkably, in the limit L — > oo the values 
Wfsa{0) and Wsa(O) coincide (within errorbars) 
which confirms the reliability of our fitting procedure. 
Results for both algorithms imply nonzero difference 
between fc and be values of the propagators, this dif- 
ference being ~ 15-^20% in the thermodynamic limit. 
Thus effect is very strong. 

Note that the asymptotic value W{0, L — oo) de- 
pends weakly (if any) on the value of the spacing a. 

In Fig. [3] we show the momentum dependence of our 
be gluon propagator D(j)) for three different values of 
/3 (i.e., for three different values of spacing a). In all 
three cases the physical volumes are approximately 
equal with aL ~ 6.0 fm. One can see that the finite- 
spacing effects are very small if to compare (3 = 7.09 
and /3 = 10.21 and are - 1 -=- 2% for |p|<400 Mev. 
For larger values of \p\ finite-spacing effects are even 
much less. We conclude that (at least) for /3 = 10.21 
we can speak about the propagator in the continuum 
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FIG. 3: be D{p) as a function of |p| for P = 4.24, (3 = 7.09 
and /3 = 10.21. In all three cases aL ~ 6.0 fm. Data for 
B = 4.24 are taken from Qll. 
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limit. 

Note that the propagator has a maximum at 
nonzero value of momentum \p\ ~ 400 MeV. There- 
fore, the behavior of D{p) in the deep infrared region 
is inconsistent with a simple pole-type dependence. 

In Fig. [4] we compare the momentum dependence 
of the be gluon propagator calculated for four differ- 
ent volumes for /3 — 10.21. Apart from p = case, 
the finite- volume dependence can be seen for compar- 
atively small values of momenta, i.e., |p|^0.5 GeV. 
For larger values of momenta the volume dependence 
quickly disappears. 

To compare Gribov copy effects for different values 



of L and various momenta we define the Gribov copy 
sensitivity parameter A(p) = A(p; L) as a normalized 
difference of the fc and be gluon propagators 



A(p) 



Df%p) - D'"'{p) 



(14) 



where the numerator is the average of the differences 
between fc and be propagators calculated for every 
configuration and normalized with the be (averaged) 
propagator. 

In Fig. \5\ we show the momentum dependence of 
A(p) for f3 = 7.09 and [3 = 10.21 for various vol- 
umes. As one can see, the Gribov copy influence is 
very strong in deep infrared and depends weakly on 
lattice spacing a. For a given value of L, the param- 
eter A{p) decreases quickly with an increase of the 
momentum. It is important that for a fixed nonzero 
physical momentum A(p) tends to decrease with in- 
creasing L. Both observations are in agreement with 
the observations made earlier in [18| and for the four- 
dimensional SU{2) theory [3l| . 

We should emphasize that our results for p = Q 
for all lattice spacings (see Fig. [T] and Fig. [5]) as well 
as our results for small nonzero physical momentum 
presented in [18[ imply that there are comparatively 
small nonzero momenta where Gribov copy effects also 
survive in the thermodynamical limit. 

The last observation is essential also for the calcula- 
tion of, e.g., screening masses in Ad theory at nonzero 
temperature where the momentum dependence of the 
gluon propagator D[p) in the infrared region is im- 
portant. 



IV. CONCLUSIONS 

We investigated numerically the Landau gauge 
gluon propagator D{p) in the 2>d pure gauge SU{2) 
lattice theory. We have employed lattices with dif- 
ferent values of L for /3 = 7.09 (a = 0.094 fm) and 
/3 = 10.21 (a = 0.063 fm). This work is the contin- 
uation of our previous paper [l^ where the most of 
calculations has been done for /3 = 4.24 (a = 0.168 
fm). 

The main goal of this work was to confirm our ob- 
servations made earlier in 18] employing larger values 
of (3 and to draw some definite conclusions about the 
continuum limit of the theory. 

The special attention in this study has been paid 
to the dependence on the choice of Gribov copies. To 
this purpose we have generated up to 280 gauge copies 
for every configuration. Our be FSA method provides 
systematically higher values of the gauge fixing func- 
tional as compared to the fe FSA and be SA methods. 
We stress that the choice of the efficient gauge fixing 
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FIG. 5: A(p) as a function of p for /? = 7.09 (Left) and for /3 = 10.21 (Right). 



procedure is of crucial importance in the study of the 
gluon propagator in the Landau gauge. 

Our main results are the following. 

1. In the thermodynamic limit L — )■ oo the value of 
D(0; L) differs from zero. This is in agreement 
with RGZ-scenario and the decoupling solution 
of DSE, and confirms the results of numerical 
computations in Ref . [ij, [l^ . 

With decreasing the lattice spacing a the value 
of D{Q) does not show the tendency to decrease. 

2. The Gribov copy effects are very strong in the 
deep infrared region. Moreover, /c propagators 
do not coincide with the be propagators even in 
the infinite- volume limit L — >■ oo (with difference 
up to- 15^20%). 

Therefore, the expectation values over the Gri- 
bov region 51 are different in the infrared from 
that calculated in the fundamental modular re- 



gion, I.e. 



(0)n ^ (0)r 



(15) 



that does not confirm the statements made in 

The finite-spacing effects appear to be rather 
small (if to compare /3 = 7.09 and /3 = 10.21). 
In the deep infrared the Gribov copy effects for 
D{p) depend weakly (if any) on the lattice spac- 
ing. So, we conclude that the difference between 
averaging over Gribov region J7 and fundamental 
modular region T persists also in the continuum 
limit. 
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